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Charge-density-wave formation in a half-filled fermion-boson transport model: 

a projective renormalization approach 
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We study the metal-insulator transition in a very general two-channel transport model, where 
charge carriers are coupled to a correlated background medium. The fluctuations of the back- 
ground were described as bosonic excitations, having the ability to relax. Employing an analytical 
projector-based renormalization technique, we calculate the ground-state and spectral properties of 
this fermion-boson model and corroborate recent numerical results, which indicate — in dependence 
on the 'stiffness' of the background medium — a Luttinger-liquid to charge- density- wave transition 
for the one-dimensional half-filled band case. In particular, we determine the renormalized electron 
and boson dispersion relations and show that the quantum phase transition is not triggered by a 
softening of the boson modes. Thus the charge density wave is different in nature from an usual 
Peierls distorted state. 

PACS numbers: 71.10.-w,71.30.+h,71.10.Fd,71.10.Hf 



I. INTRODUCTION 

Charge density waves (CDWs) are broken symmetry 
states of metals that predominantly appear in materi- 
als which have a highly anisotropic crystal and electronic 
structure^. The formation of CDWs strongly depends on 
the band-filling and on the topology of the Fermi sur- 
face. Concerning the latter, one-dimensional (ID) sys- 
tems are peculiar, because their Fermi surface consists 
of two points only. Not surprisingly, electron-electron 
and electron-phonon interactions, which are the driv- 
ing forces behind most metal-insulator transitions, have 
more impact in reduced dimensions. At the same time, 
however, quantum fluctuations and finite-temperature 
effects — both counteracting any development of long- 
range order — become increasingly important as well. 
Then, all in all, in ID, usually a rather complex interplay 
between the charge, spin, orbital and lattice degrees of 
freedom evolves, which, of course, will strongly affect the 
transport properties of the system. Prominent examples 
are quasi- ID halogen-bridged transition metal complexes 
(MX-chains)2. 

The theoretical description of such highly correlated 
ID systems is frequently based on 'microscopic' mod- 
els like the half-filled SSH 3 , HolsteiniA Peierls!, Hub- 
bard£, (quarter- filled) t-J models 9 , or combination of 
thesei£. Thereby, the complexity of the electron-electron 
or electron-phonon interactions in the Hamiltonians usu- 
ally prevents the exact (numerical) solution of the model 
in the thermodynamic limit, which would be necessary 
in order to pinpoint a true quantum phase transition be- 
tween a metal and a CDW. 

To some extent, the state of affairs improves if one con- 
siders simplified transport models instead, where particle 
motion takes place in an effective background medium. 
The 'background' reflects the correlations inherent in the 



system, e.g., the charge-, orbital- or spin-order in a solid. 
That quasiparticles move through an ordered insulator is 
a very general situation in condensed matter physic o 11 ' 12 . 
This scenario also applies to soft matter systems like 
DNA, where the charge transport on the backbone is 
affected by the 'configuration' of chemical side groups, 
which, vice versa, depends on the physical presence of 
charge carriers^ 3 -. 

Along this line, a novel fermion-boson quantum trans- 
port model, 
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has been proposed a few years ago^, and was shortly 
afterwards solved for a single particle (N e — 1) in a ID 
infinite system^ by a variational numerical diagonaliza- 
tion technique^. The Hamiltonian (fTJ) mimics the 'back- 
ground' by bosonic degrees of freedom [b\ ], which in- 
fluence and even may control the transport of fermionic 
particles [c,-^] on the N sites of a regular lattice. Every 
time a particle hops between nearest-neighbor Wannier 
sites it creates (or destroys) a local excitation of 

energy w& in the background medium at the site it leaves 
(it enters). Clearly these distortions tie the particle to 
its origin (cf. the string effect as a hole moves in a Neel 
spin background). Of course, any distortion of the back- 
ground can heal out by quantum fluctuations (again, one 
can have spin fluctuations in mind). Accordingly, in the 
Hamiltonian (TTJ), the Af,-term was included, allowing for 
spontaneous boson creation and annihilation processes. 

Strong correlations, nevertheless, may evolve in a sys- 
tem described by the Hamiltonian (JTJ) , provided the back- 
ground excitations have a rather large energy and the 
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ability of the background medium to relax is small, i.e., 

Q = — > 1 and A = — < 1 . (2) 

tb tb 

For the half-filled band sector (N e = N/2), these correla- 
tions may even drive the system into an insulating state 
by establishing CDW long-range order. This has been 
shown quite recently for the ID case: small cluster diag- 
onalizations— and density matrix renormalization group 
(DMRG) calculations supplemented by finite-size scal- 
in g 17 i ls give strong evidence for a Tomonaga-Luttinger- 
liquid (TLL) CDW quantum phase transition as A be- 
comes small at large enough Q. These purely numeri- 
cal approaches rely on an (inevitable) truncation of the 
bosonic Hilbcrt space. Determining the metal-insulator 
phase boundary this seems to be uncritical, because the 
CDW found at half-filling is a few-boson state. The situ- 
ation becomes more difficult if we enter the fluctuation- 
dominated regime of small tt, where many bosons are 
excited in the system. 

In the present work, we investigate the fermion-boson 
transport model (fTJ) by means of an analytical approach, 
which avoids these disadvantages. This approach, called 
projective renormalization method (PRM)i2, is based on 
a sequence of discrete unitary transformations, so that — 
in contrast to continuous (e.g. flow-equation based) uni- 
tary transformation schemes^ — a direct link to pertur- 
bation theory can be provided. The method has already 
been successively applied to a number of many-particle 
model o 19 ' 21 " — . Here we will analyze the ground-state 
and spectral properties of the Hamiltonian ((T|) exclusively 
for the half-filled band case, in both the metallic and in- 
sulating regimes. In particular we study the signatures 
of the TLL-CDW transition in terms of the renormalized 
quasiparticle band and boson dispersion, and the boson 
spectral function. The paper is organized as follows. In 
Sec. II. A we briefly resume the basic concepts behind 
the PRM approach. The application of the PRM to the 
fermion-boson transport model will be described in de- 
tail in Sec. II B. Section III presents the results of the 
numerical evaluation of the renormalization equations. 
We conclude in Sec. VI. 



II. THEORETICAL APPROACH 

A. Projector-based renormalization method 

The PRM starts from the usual decomposition of a 
many-particle Hamiltonian into a solvable unperturbed 
part Hq and a perturbation Hi, where Hi should not 
contain any part that commutes with Hq. Thus, the per- 
turbation Hi consists of transitions between the eigen- 
states of Hq with non-vanishing transition energies. The 
basic idea of the PRM is to construct an effective Hamil- 
tonian H\ = Hq x + Hi \ with renormalized parts -£/o,a 
and Hi, a, where all transitions with energies \E£ x — E™ x \ 



larger than a given cutoff energy A are eliminated. Efi x 
and E™ x denote the eigenenergies of i?o,A- 

The renormalization procedure starts from the cutoff 
energy A = A of the original model H and proceeds in 
steps of AA to lower values of A. Every renormalization 
step is performed by means of a unitary transformation, 

H x _ AX =e x ^H x e~ x ^\ (3) 

The generator Xa,aa of the unitary transformation has 
to be fixed appropriately (for details see Ref.Ql). For 
instance, in lowest order perturbation theory, it reads 

^a,aa = t Qa-aa-Hi.a • (4) 

-L^A 

Here, Lo,a is the Liouville superoperator of the 'un- 
perturbed' Hamiltonian H \ , which is defined by the 
commutator of -ffo.A with any operator variable A, i.e. 
Lo,a^4 = [Hq^x, A], and Qa-aa is a projection superoper- 
ator, which projects on all transitions with respect to the 
eigenspectrum of i?o : A-AA with transition energies larger 
than A — AA. In this way difference equations can be de- 
rived which connect the parameters of H\ with those of 
H\-a\i and which are called renormalization equations. 

The limit A — >• provides the desired effective Hamilto- 
nian H = H\^,q = Hq.x^q where the elimination of the 
transitions originating from the perturbation Hi leads to 
a renormalization of the parameters of H. Note that H 
is diagonal or at least quasi-diagonal and allows to eval- 
uate physical quantities. The final results depend on the 
parameter values of the original Hamiltonian H. Finally, 
we note that H and H have the same eigenvalue prob- 
lem since both Hamiltonians are connected by a unitary 
transformation. 

To evaluate expectation values of operators A, formed 
with the full Hamiltonian, we have to apply the unitary 
transformation as well, 

TiAer m 

where we define A{\) = e Xx Ae~ Xx and A = A(X -> 0). 
Thus additional renormalization equations are required 
for A(\). 



B. Application to the two-channel transport model 

1. Renormalization equations 

We first rewrite the model ([TJ, performing a unitary 
transformation bi i— > bi + Xb /u)b that eliminates the bo- 
son relaxation term in favor of a free-particle hopping 
channel, 
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with 



11 = 2^ = 2^ 



(7) 



This makes the two transport channels contained in the 
Hamiltonian ([1} explicit: the coherent particle transfer, 
which takes place on an energy scale (oc i/), and the 
boson- affected hopping (oc 

Next, in order to exploit the translation invariance, we 
consider the Hamiltonian ([5]) in momentum space 



k q 



^Y,9k(bl c l c k+ q + b q c l+ q c k) 

k,q 



(8) 



In what follows, we consider a ID lattice with lattice 
constant a, i.e., e k = —2tfCOska and g k = —2t b coska. 

Going forward, it turns out to be useful to remove the 
mean-field part from the fermion-boson coupling term. 
Defining fluctuation operators, 

5(c{ C k+q) = c l C k+q ~ ( c l C k) tqfi , (9) 

the Hamiltonian © takes the form 

H = J2 £ k c i c k+ uj bY. b l b i 

k q 

+ -7^Y,9Mc k )( bt o + b o) 

V k 

+ -7^ E 9k W{c{c k+q ) + b q S(cl +q c k )] . (10) 

k,q 

Obviously, the purely bosonic part of H, i.e. the second 
and third term of Eq. (ITU)) , can be diagonalized by a shift 
of the bosonic operators. Introducing new bosonic cre- 
ation operators, 



B l= b l 



N v Ub 



(ii) 



the Hamiltonian (fTU|) can be rewritten as H = Hq + Hi 
with 



T. B l B q + ^-{H9k{c\c k ))\ (12) 



R i = -J^T,9k[Bl6(clc k+q ) + B q 5{cl +q c k )].(Vi) 

k,q 

Following the ideas of the PRM approach, we make the 
following ansatz for the renormalized Hamiltonian H\ 
(after all transitions with energies larger than A have been 



integrated out), H\ — H Q ,\ + Hi t \, where 

H ,x = J2e k , x clc k + J2 A k,xcl c k+Q 

k k 

+ y ^u qtX B\B q + E x , 
q 

H ^ = ^T,9kQk^)[BlS(clc k+q ) + B q S(cl +q 

k,q 

The 9-function in SB 



(14) 



3 C fe) 



(15) 



6fc, 9 (A) = 6(A - |£fe,A - £ k +q,\ + Wq,\\) 

guarantees that only transitions with excitation energies 
smaller than A remain in H\ \. In H$\, also a symmetry 
breaking field, A kt \, was introduced which couples to- 
gether particle-hole excitations with wave vectors k and 
k + Q, where Q = ir/a. Note that the renormalization of 
Afe.A may lead to a transition to a CDW ground state at 
half filling. 

By integrating out all transitions between the cutoff 
of the original model A and A = 0, all parameters of the 
original model will become renormalized. To find their A- 
dependence, we derive renormalization equations for the 
parameters s k> \, A kj \, w q ,\, and E\. The coupling pa- 
rameter g k is not renormalized when we restrict ourselves 
to lowest order perturbation theory in one renormaliza- 
tion step. The initial parameter values are determined 
by the original model (A = A): 



-fc,A 



s k -2g k -J2^(cU k >) 

k' 



1 / 



y^9k{c\c k ) 



A i. 



k.X 



, and uj t 



Ub . 



(16) 

(17) 
(18) 



Let us assume that the symmetry breaking field ~ A^a 
can be considered as small compared to the hopping part 
in Hq \. In this case, the dynamics of -£/o,a is approxi- 
mately governed by 



[H ,x,cl] = £ k ,\c\, 



(19) 
(20) 



Following [19f . the lowest order expression generator 
X\ aa is obtained as 



X 



—y 

y/N kq £fc,A - £fe+g,A + ^q,\ 



A,AA 



g k & k , q (\,A\) 



[B\6{clc k+q ) - B q S(cl +q c k )] . (21) 



Here, 



e k , q (X, AA) = 6(A - |e fe , A - £ k+q ,x + u q<x \) (22) 
x6(|e fe ,A-AA - £ k +q,\-Ax + w 9 ,a-aa| - (A - AA)) 
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is a product of two 0-functions which assure that only 
excitations between A and A — AA are eliminated by the 
unitary transformation (|3]). From Eq.Q, the Hamilto- 
nian -Ha- a a is easily evaluated within second order per- 
turbation theory, 



H\-a\ — H\ + [X\.ax, Ho,\] 



(23) 



An alternative expression for Hx-ax is obtained by 
replacing A in Eqs. (fl"4|) . (fT5|) by the reduced cutoff A— A A. 

Thus, i?A-AA = #0,A-AA + #1,A-AA with 

H ,X-AX = J2 £k > X - AxC k C k+J2 A k,\-A\clc k+Q 



H\ X—AX — 



1 

-L^ 5fc e M (A-AA) 
x[B\5{c\c k+q ) + B q 5{cl +q c k ) 



(24) 



(25) 



A comparison of expression ([23]) with ([24]) , (|25[) leads to 
the renormalization equations which connect the param- 
eter values of the Hamiltonian at cutoff A with those at 
cutoff A — AA. We obtain 



l sri(B^ c x gfcQ M (A,AA) fl c g ^_ g 8 fc - g , g (A,AA) 

£fe,A-AA - £fc,A = T7 X I < + " T + < ~ + 1 " ) > 26 

1 W c c \ gi®k,q{\ AA) 
W g ,A-AA-W g ,A - -jy 2^("fc ~ "fc+g h XT—' ( 2? ) 

a a 1 / g 2 k e k>q {\)G k+Q>q (\,A\) g 2 k+q Q k+q ,- q (X)e k+q+Q ^ q (X,AX) 

A k< X-AX - A fcj A = — > < ; 1 ; 

g I £k +Q> X ~ £ k+q+Q,X + ^q,X £k+q+Q,X ~ £ k+Q,X + ^>q.X 



g g9 fc+q , g (A)e fc , g (A,AA) | aj +q Q k+q+Q ,- q [\)Q k+q ,- q {\, AA) 

£fe,A — £fe+g,A + ^ g ,A £ k+q,X — £k,X + ^>q,X 

1 /e fc+Q>Q (A)e fc ,, Q (A,AA) e fc . Q (A)e fc , + Q, Q (A,AA) 



a k+q 



^ y I £ k',A^£)c'+(3,A+ w Q,A £fc'+Q,A — £fe',A + ^Q,A 

9 fc ,, Q (A)9 fc , Q (A,AA) e,> +Q ,Q(A)9 fc+ Q, Q (A,AA) | 
; ; f 9k9k'd k , 

£k,X — £k+Q,X + WQ,A £fe+Q,A — £/c,A + ^q,a J 



(28) 



Since the equation for the energy shift Ex- ax is not 
needed in the following it has been left out for brief- 
ness. In Eqs. ([26| - ([28|) . we have defined new expectation 
values n% = (c{c k ),nf = (B^B q ), and d% = (c{c k+Q ), 
which arc formed with the full Hamiltonian H. Suppose 
these expectation values are known, the renormalization 
between the cutoff A of the original Hamiltonian H and 
A = leads to the Hamiltonian H, 

H = J2 tkcW+y A k c{c k+Q +y u q B\B q +E , (29) 

k k q 



where i k , A k , u q and E denote the parameter values at 
A = 0. 

Note that the fully renormalized Hamiltonian H de- 
scribes an uncoupled system of renormalized (dressed) 
electrons and bosons. Both parts are quadratic either in 
the fermionic or in bosonic operators. By a rotation in 
the fermionic subspace the electronic part of H can eas- 
ily be diagonalized. Thus, any expectation value can be 
evaluated. This property will be used in the following in 
order to evaluate the yet unknown expectation values n k , 
n q , and d%. 
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2. Expectation values 



and 



The expectation values can be evaluated self- 
consistently within the PRM formalism by applying the 
same unitary transformation as was used before for the 
Hamiltonian. Following Eq.Q, for instance, n k can 
be expressed by n k = (c£(A)c*(A))h a , where {■■■)h x 
means the average formed with H\ and c£(A) is given 
by <4(A) = e Xx c\e~ Xx . For the transformed operators 
<4(A) and Bj(A) we use the ansatz 

4(A) = a k , X c\ + \h, q ,\ B q A+ q + lk,q,\ B t C l-q 



(30) 



flj (A) = 4> q X Bl + %tX B- q + J2 4>k, q ^(cl +q c k ) , (31) 



respectively. The operator structure is again taken over 
from the lowest order expansion of the unitary transfor- 
mation, apart from the second term in -Bj(A) which is due 
to higher-order terms. For the A-dependent coefficients 
otk,X, Pk,q,X, • • • also renormalization equations have to be 
derived, 



Oik, X— AX — OLk,\ 



2^ 



9k 



Pk,q,X-AX — Pk,q,X = 

lk, q ,X-AX — lk,q,X — 

4>q,X-AX — <j)q.X = 

Vq.X-AX - %A = 

i>k,q,\-AX — ~>Pk,q,X — 



2N 

1 



£k,X — £k+q,X + LOq.X 
9k-q 



a fe , A e M (A,AA) 



9k 



'N £k,x ~ £k+q,X + Wg,A 
1 9k+q 



£k-q,X — £k,X + ^q,X 

a.k,x®k, q {\ AA) , 



afe,A@fc-g,g(A, AA) , 



N £k+q,X — £k,X + <^q,X 



a*: : A@fe+ g ,-g(A, AA) , 



2iV 



N £-~t 



l k+q) 



9k 



£k,X — £k+q,X + Mq,\ 



<i>q,X®k,q{\ AA) , 



27V 



1 



.9/v 



£k+q,X ~ £k,X + Wg,A 



%,A©fe+9 ,-g(A, AA) , 



iV l£k,X~£k+q,X+Ulq,X 



0q,X®k,q( X , AX ) - 



9k+q 



£k+q,X~ £fc,A+Wq,A 



Vq,xQk+q ,-g(A, AA) 



(32) 

(33) 
(34) 

(35) 

(36) 

.(37) 



Integrating these equations between A (where a k A = 
<p q ^ = l and all other coefficients zero) and A = 0, we 
arrive at the final result for n k , n B , and d k : 



Hi = ain 



k IL k 



[Kqi 1 + "f )"*+g + llq^qnl-q 



B T2~B 

n q = 4> q n q 



+fj 2 q (l + n B q ) + J2^lqn c k+ q^ ~ n%), 

k 



d k = Oka k +Qd c k 



lk, q lk+Q, q n q 



(38) 



(39) 



(40) 



l k—q ■ 



Here, a^, Pk.q, ■ ■ • denote the fully renormalized parame- 



ter values at A = 0. Similarly, n k ,n B and d k are expec- 



tation values defined with H, i.e., 



n k = { c k c k)f[, 
^fc = ( c I C fc+Q)_H- ' 



(41) 
(42) 
(43) 



where the fully renormalized Hamiltonian H is given by 
Eq.©. 



5. Dynamical correlation functions 
Let us consider the boson spectral function, 

W = 5— / <[&*(*), e^dt, (44) 



6 



and the two electronic one-particle spectral functions 



i 

2^ 



c fc (i) c\ 



1 /"°° / 



dt. 



dt. 



(45) 
(46) 



Here, A k (u>) describes the creation of an electron with 
momentum k at time zero and its annihilation at time 
t whereas in A k (uj) first an electron is annihilated. As 

I 



it is well-known, A k (oj) and A k (oj) can be measured by 
inverse photoemission and by photoemission. 



To evaluate Eqs. (|44J)— (J46J) within the PRM approach, 
we use again that expectation values are invariant with 
respect to a unitary transformation under the trace. 
Thus, C q (oj), A k (ui), and A k (oj) can easily be computed 
if the bosonic and electronic one-particle operators are 
transformed in the same way as the Hamiltonian. In this 
way we obtain 



; 2 ~2 

= ^d( w -a,) + ^( w + o_,) + £$„^ 



^-d(e k+q -£k- w), 



£k^ 



CM 

= a 2 k n c k S(uj - s k ) [Pk, q ( l + nf) n c k+q S(uj + Q q - e k+q ) + % q nf n%_ q 8{ui - ui q - £ fe _ 9 ) 

A+(w) = &l{l-nl)5{u-s k ) 



+ 



(47) 



(48) 



(49) 



where terms with two bosonic creation or annihilation 
operators have been neglected. The 4> q , fj q , and tjjk,q 
are the zero-A coefficients taken from the evaluation of 
Eqs.([25l)-(I371). 

Let us emphasize that the expressions (|44)) . (|45l and 
(|ini fulfill the sum rules, 

/oo 
dojLoC q (oj) = 1 (50) 

and 

/oo 
doj[A+(u J )+A-(o J )] = l, (51) 
-OO 

respectively, which also hold if Eq. (|4"T1) for C 9 (w) and 
the corresponding expressions for A k (oo) and ^4^(w) are 
inserted. 



4- Numerical analysis 

The set of renormalization equations (|26[) — (f28]) and 
(I5^1) - ([57|) has to be solved numerically. To this end, we 
choose some initial values for the expectation values en- 
tering the renormalization equations. Using this set of 
quantities, the numerical evaluation starts from the cut- 
off A of the original model H and proceeds step by step 
to A = 0. For this procedure we consider a lattice of 
N = 500 sites in one dimension. The width of the en- 
ergy shell AA was taken to be somewhat smaller than the 
typical smallest energy spacing of the eigenstates of i?o,A- 
For A = 0, the Hamiltonian and the one-particle opera- 
tors are fully renormalized. The case A = allows the 
re-calculation of all expectation values, and the renor- 
malization procedure starts again with the improved ex- 
pectation values by reducing again the cutoff from A to 
A = 0. After a sufficient number of such cycles, the ex- 
pectation values are converged and the renormalization 
equations are solved self-consistently. Convergence is as- 
sumed to be achieved if all quantities are determined with 
a relative error less than 10 -5 . The dynamical correla- 
tion functions (|4T)) -(|4*^ | are evaluated using a broadening 
in energy space that is equal to AA. 



III. GROUND-STATE PROPERTIES 



A. DMRG phase diagram 



In order to classify the PRM ground-state and spec- 
tral properties given below, we first present in Fig. Q] 
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a refined version of the DMRG ground-state phase di- 
agram of the half- filled fermion-boson model ([1}. Here 
the phase boundary, separating the insulating phase with 
CDW long-range order from the metallic TLL phase in 
the A-f2 plane, was obtained from the N — > oo extrap- 
olated values of the Luttinger liquid parameter K p and 
the single particle (charge) gap A ^ 17 ! 24 . In the limit of 
large CI, the background fluctuations, associated with any 
particle hop, are energetically costly. As a result the mo- 
tion of the particle is hindered and charge ordering be- 
comes favorable if A, describing the ability of the back- 
ground to relax, is sufficiently low. By contrast, for large 
A (A > A c (Cl — oo) ~ 0.1588), we find metallic behavior 
for all CI. In the limit of small CI, the rate of bosonic fluc- 
tuations (oc Cl^ 1 ) is high. Then, in no way, correlations 
emerge within the background medium. The DMRG re- 
sults suggest that for A = 0, i.e. when the relaxation 
channel is closed, the ground state is nevertheless metal- 
lic below a finite critical boson energy Cl c (A = 0). Let 
us re-emphasize that coherent particle hopping is possi- 
ble even when A = 0, due to a six-step vacuum-restoring 
hopping process*^, 



R 



(6) 

i,i+2 



^1+2^1+1^1^+2^+1-^ 



(52) 



J, 



with R\ = c\c i+1 b l and h\ = cjq^. r\% 2 
'effective' (coherent) next nearest-neighbor transfer. 
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FIG. 1: (Color online) DMRG phase diagram of the two- 
channel fermion-boson transport model JTJ [see also Eq. ([6|] 
for the ID half-filled band case (the dot-dashed line is a guide 
to the eye). The crosses mark the parameter values used 
within our PRM calculation when passing the TLL— >CDW 
transition at fixed O = 10 {a} and A = 0.001 {b} (green 
arrows), respectively. 

Concerning the nature of the metal-insulator quan- 
tum phase transition, it is a moot point, whether the 
TLL-CDW crossover in the half-filled fermion-boson 
model (HJ), taking place at relatively large Cl, bears some 
resemblance to the usual Peierls transition in the spinless 
fermion Holstein model^ i 18 ' 22 . In this regard the ques- 
tion of boson softening will certainly be of importance. 



B. CDW order parameter 

To analyze the nature of the metal-insulator transition 
of ([1]) in more detail, we calculate in the following a set 
of characteristic quantities by the PRM for the (ID half- 
filled) infinite system. Thereby we cross the TLL— ^CDW 
transition in the following figures at fixed CI [panels (a)] 
and A [panels (b)] (cf. Fig. Q] lines {a} and {b}, respec- 
tively). As for the half-filled Holstein model 2 !, the CDW 
structure of the insulating state shows up in the correla- 
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FIG. 2: (Color online) Zero-temperature expectation value 
d%, indicating CDW order in the half- filled fermion-boson 
model Q, at n = 10 [upper panel (a)] and A = 0.001 [lower 
panel (b)]. 

Figure [2] displays the variation of this expectation value 
when the wave vector k runs through the half ID Bril- 
louin zone. Obviously we have d\ = in the metallic 
phase (blue dot-dashed line). Entering the CDW state 
d c k acquires finite values, whereby the maximum of d\ is 
at k = ir/2 (for this case Q = ir connects both Fermi 
momenta kp — ±ir/2). If the charge order is perfect (the 
particles are localized in an A-B structure without any 
charge fluctuations, i.e. the lower and upper bands are 
flat), we find g?£ = 1/2 for all k. This tendency becomes 
apparent by comparing the results obtained in the CDW 
phase for different A (cf. upper panel A = 1/8, 1/12 
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and lower panel A = 1/1000; recall that A describes the 
ability of the background to relax.) 



C. Fermion dispersion and quasiparticle weight 

Next we investigate the renormalization of the 
fcrmionic band structure, 



E k = 



£k + £k+Q 



± 



£k — £fe+Q 



(53) 



see Fig. [3] In the metallic regime, of course, there is 
no gap at the Fermi energy (Fermi vector kp = 7r/2), 
and Afc, given in the inset, is zero for all k. While 
for 2A f2 the free transport channel (oc tf) dom- 
inates even when f2 is large (see Eq.([7])), the bosonic 
degrees of freedom will strongly affect the transport for 
small A. As a consequence, 'coherent' transport takes 
place on a strongly reduced energy scale only (we have 
tf/U = 1/30 [tf/tb = 1/400] for the blue dot-dashed line 
in panel (a) [(b)]). 






1/12 
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FIG. 3: (Color online) Dispersion of the occupied lower 
(empty upper) quasiparticle band Ek with k < it/ 2 (k > 7r/2) 
for fi = 10 [panel (a)] and A = 0.001 [panel (b)] (the Fermi 
energy sets the zero-point of energy). Note that band energies 
are differently scaled by 2tf and 2004/ in (a) and (b). The 
insets show the fc-dependence of A fe . 
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FIG. 4: (Color online) Spectral weight, a\, of the coher- 
ent (quasiparticle) part of the electronic one-particle spectral 
function (|fg|). We have Q. = 10 in panel (a) and A = 0.001 in 
panel (b). Notations are as in Fig. [3] 



The coefficient 5|, depicted in Fig. [4] gives the weight 
of the corresponding coherent part of the single-particle 
spectral function (|48|) . A^((j) can by probed by angle- 
resolved photoemission experiments. At very large f2 
(and small A), the particles will solely move by the above 
mentioned six-step process ()52j) . Then the resulting 
'quasiparticle weight', a|, is nearly one [see Fig. 2] (a)], 
and shows a very weak A-dependence. For small Q < O c 
we enter the fluctuation-dominated regime and the na- 
ture of the metallic state changes noticeably. In accor- 
dance with recent dynamical DMRG data for the single- 
particle spectra, which show that the absorption spec- 
trum is over-damped near k — 0, n because of intersect- 
ing bosonic excitations, we find af ~ 1 in the vicinity 
of k F only [cf. Fig. 3 in Ref. [U and Fig. H (b), blue 
dot-dashed line]. 

In the insulating regime, the renormalized band struc- 
ture Ek is gaped (see Fig.[3J dashed and solid lines). The 
inset in Fig. [3] clearly shows an increase of the gap at 
k = k F as A gets smaller. Note that the size of the gap 
is equal to 2Ak F ■ While Afe is symmetric around k = n, 
Ek is not. The reason is that doping a perfect CDW, 
states with one particle removed are connected by the 
six-step hopping process (j5"2"j) . whereas a two-step hop- 
ping process relates states with an additional particle^. 
In this way the collective particle-boson dynamics leads 



9 



to a more pronounced flattening of the coherent band for 
k < hp, i-e., the widths of the highest photoemission and 
lowest inverse photoemission band diffe r 16 ! 18 . It is en- 
couraging that our analytic PRM approach reproduces 
this non-trivial correlation-induced (mass-) asymmetry. 
Let us emphasize that the aj* given in Fig.|4]for the CDW 
case (dashed and solid curves) belong to the highest pho- 
toemission band in the whole interval [0,tt] (the corre- 
sponding Ek is not depicted in the region n/2 < k < tt 
in Fig. [3]). Compared to the metallic phase the spectral 
weight of the lower CDW band is significantly changed 
for intermediate-to-small boson frequencies only. 



D. Boson dispersion and occupation numbers 

The Einstein bosons, describing excitations of the 
background, gain a dispersion owing to the coupling to 
the fermions. The renormalization of the boson disper- 
sion, ujq/ujb, is displayed in Fig. [5] It is rather weak for 
large SI in both the metallic and insulating states. For 
smaller boson frequencies, we find a strong renormaliza- 
tion in the TLL phase (up to 50% for £1 = 0.8) at larger 
momenta [see dot-dashed curve in panel (b)]. This is in 
accordance with the over-damped single-particle excita- 
tions observed in the ARPES spectra^ and, of course, 
also shows up in the depletion of ot\ away from k = tt/2 
[see Fig. 01(b)]. 

Most notably, for the CDW state, we observe a hard- 
ening of the boson modes near the Brillouin zone bound- 
ary. This holds in the whole Q region and means that 
the TLL— >CDW transition is unlike the usual displacive 
Peierls transition which, in general, is accompanied by 
the softening of the q = tt boson (phonon)&22. In our 
case, the CDW state is driven by the stiffness of the back- 
ground, being most pronounced at large fl and small A. 
By contrast, when f2 is small, i.e. the background read- 
ily fluctuates, the kinetic energy part will naturally over- 
compensate any potential energy gain by charge order- 
ing. Another reason for the absence of boson softening 
might be the particular form of the fermion-boson in- 
teraction. As can be seen from the Fourier transformed 
Hamiltonian (jSJ) , the fermion-boson coupling vanishes for 
k = ±7r/2, i.e., precisely for the Fermi momenta of the 
half-filled band case. 

It may be worthwhile to demonstrate that our PRM 
approach has the advantage that all features of the results 
for Q q and all other renormalized quantities can easily 
be understood on the basis of the former renormaliza- 
tion equations. For simplicity we shall restrict ourselves 
to the renormalization of Co q in the case of large ft. In 
this regime, from Eq. (|27p one may point out the stiffen- 
ing of the boson modes. Since the boson energy ujb is 
much larger than the electronic bandwidth, for all A a 
positive energy denominator (ek,\ — £k+q,\ + ^q,\) i s ob- 
tained. Nevertheless, in the k sum on the right hand side 
of Eq. (|27p there are as many negative as positive terms 



1.002 





1/12 


A = 


1/8 


A = 


1/6 



1.001 



3 

a 

;3 



0.999, 



1.2 




0.8 



'3 



0.6 



0.4, 



(b) 



— a = 2.0 

■-■ £2 = 1.8 
-■- Q = 0.8 



0.25 



0.5 

q / n 



0.75 



due to the factor [n\ 



L k+q 



). Since from (n c k 



l k+q 



)<0 



FIG. 5: (Color online) Renormalized boson frequencies ui q for 
the fermion-boson transport model (flj with fi = 10 [panel 
(a)] and A = 0.001 [panel (b)]. 



it follows that {£k,\—£k+q,\) > 0, the negative terms have 
larger energy denominators and are always smaller than 
the positive terms. The resulting renormalization of uj q is 
therefore positive for all q values and largest for q — tt due 
to the smallest energy denominator. Furthermore, since 
gk oc tb and uib is large the renormalization contributions 
in Eq. (|2"T|) are of the order of t\l<^b = ^b/fl 2 <C oJb which 
gives rise to the weak dispersion of uj q observed in Fig. [5] 
For smaller values of f2 the bosonic and fermionic energy 
values in the denominator of Eq. f|27fl can become com- 
parable which immediately leads to a strong dispersion 
of ui q (see dot-dashed curve in Fig. [5] and solid curve in 
Fig. [5J . 

Figure [6] gives the (g-resolved) boson occupation num- 
bers. As one can see from Eq. ([59"]) . this quantity 
for T — acquires finite values solely by coupling to 
fermionic degrees of freedom. Note that the first term in 
Eq. (|39p vanishes for T — 0. We see that the formation 
of the CDW state is accompanied by a finite occupation 
value of the q = boson mode, which is about two orders 
of magnitude larger if one compares Uq for the CDWs 
established at fl = 1.8 and f2 = 10, respectively. Dif- 
ferent from the Holstein-model CDW (Peierls) phased, 
the CDW phase of the half-filled fermion-boson transport 
model ([T]) is always a few-boson state however. Referring 
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FIG. 6: (Color online) Transformed boson expectation value 
nf = (B\B q ) at n = 10 (a) and A = 0.001 (b). 



to this, our PRM results corroborate previous small clus- 
ter exact diagonalization data^. As can be seen from the 
third term of Eq. (1391) . bosons having finite momentum 
give rise to an effective fermion interaction on neighbor- 
ing sites. 



IV. SPECTRAL PROPERTIES 

While the fermionic single-particle spectral function of 
the transport model ([T]) was previously calculated for fi- 
nite clusters by exact diagonalization 1 ^ and dynamical 
DMRG 1 ® techniques, the spectral response of the bosons 
has not been studied so far. The PRM allows to inves- 
tigate the interrelation between fermion and boson dy- 
namics by computing the boson spectral function, C q (ui), 
according to Eq. (|4"T)l , for the ID infinite system^. 

Figure [7] shows C q (uj) in the metallic regime, for dif- 
ferent fi and A parameters. For very large f2 = 10 
[panel (a)], the boson energy is hardly renormalized by 
the coupling to the fermions. Accordingly we observe a 
strong signal at the bare boson frequency fi/2i/ = 100 
(first term in Eq. (|47p: the second term in Eq. (|4"?| 
will not contribute because there are no states available 
with uj = -Q/2i/). The third term in Eq. (gTJ) detects 



particle-hole excitations and leads to the two incoherent 
absorption bands in C q (uj) running from q = — > n with 
energies between uj = and uj ~ ±10i/ = ±4/2. At 
small = 0.006, see panel (b), the (one-) boson excita- 
tion is located within the fermionic band. As a result we 
find a strong renormalization of the bare boson frequency 
(see also Co q in Fig. [5]), leading to the dispersive signal in 
the range ui/2tf ~ 0.5 x 10 -3 . . . 10~ 3 . We note that for 
A = 0.01 used in panel (b) the fermion-boson coupling 
is small in comparison with the free fermion bandwidth 
(we have gk/£k = h/tf = 0.3 in the model ©), hence 
the effect of multi-boson absorption processes is negligi- 
ble. The lower two panels of Fig. [7] demonstrate how A 
affects the boson absorption at fixed f2. In panel (c), for 
A = 1 and fl = 4, the boson frequency is larger by a fac- 
tor of two than the 'free' fermion bandwidth (4ty = 2), 
whereas they have the same size for the C q {uj) spectrum 
with A = 2 shown in panel (d). Quite differently, in the 
former case, the bare boson mode hardens, while it soft- 
ens near Q — n in the latter case, where the fermion and 
boson degrees of freedom are strongly mixed. This be- 
comes even more visible by comparing the corresponding 
(dashed and solid) ui q curves in Fig. [5] 



V. CONCLUSIONS 

To summarize, we adapted the projective renormal- 
ization method to the investigation of a novel two- 
channel fermion-boson model, describing charge trans- 
port within a background medium. By large-scale nu- 
merical DMRG studies this model has been proven to 
show a metal-insulator quantum phase transition for the 
one-dimensional half-filled band case. The transition is 
triggered by strong correlations evolving in the back- 
ground and typifies as a Luttinger-liquid charge-density- 
wave crossover. 

Our analytical approach captures this TLL-CDW tran- 
sition for the infinite system, without truncating the 
bosonic Hilbert space as in purely numerical investiga- 
tion schemes. Therefore, the PRM is particularly well 
suited to analyze the bosonic degrees of freedom when 
passing the metal-insulator phase boundary. 

In the course of the renormalization procedure of the 
fermion-boson Hamiltonian we end up with a model of 
noninteracting — but dressed — electrons and bosons. In 
this way, the renormalization of the fermion band disper- 
sion and of the boson frequency is obtained, and various 
ground-state expectation values were calculated. More- 
over, we derived analytical expressions for the single- 
particle (inverse) photoemission spectra and for the bo- 
son spectral function, which allows us to pinpoint the 
most important absorption and emission processes during 
particle transport, and throws some light on the nature 
of the TLL-CDW transition. 

In particular, we show that the insulating CDW phase, 
realized for large boson frequencies Q and small boson re- 
laxation parameter A, is characterized by a gapful, mass 
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FIG. 7: Boson spectral function C q (uS) of the half-filled fermion-boson transport model {TJ with: SI = 10, A = 1/4 [panel (a)]; 
fi = 1/150, A = 1/90 [panel (b)]; f2 = 4, A = 1 [panel (c)]; Q, = 4, A = 2 [panel (d)]. The frequency tv is given in units of 2tf. 
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FIG. 8: (Color online) Renormalized boson frequencies for the 
parameters used in Fig. [7] CI = 10, A = 1/4 (green double- 
dot-dashed line; CI = 1/150, A = 1/90 (red dot-dashed line); 
CI = 4, A = 1 (blue dashed line), CI = 4, A = 2 (black solid 
line) . 



by numerical approaches. Here we observe a well-defined 
electron band in the vicinity of kp only. Due to many 
intersecting boson branches the (inverse) photoemission 
spectra are 'overdamped' near the Brillouin zone bound- 
aries. 

The boson spectral function and renormalized boson 
dispersion clearly indicate that the TLL-CDW transition 
is not accompanied by a softening of the zone-boundary 
boson mode. Rather a harding of the Q = it boson is ob- 
served. This might be partially attributed to the vanish- 
ing Fourier-transformed fermion-boson coupling term at 
wave-numbers ±7r/2, which denote the two Fermi points 
for the half- filled band case. The situation changes if we 
look for a metal-insulator transition for other commensu- 
rate band-filling factors, e.g., at quarter filling. Whether 
the system there undergoes a soft-mode transition for 
small CI is an interesting open question that deserves fu- 
ture efforts. In this connection other kinds of charge- 
ordered states should be also considered. 



asymmetric band structure. Thereby the lower occupied 
band is almost fiat by reason that transport is only pos- 
sible through a vacuum restoring six-step hopping pro- 
cess. The CDW phase is a few-boson state. By contrast, 
the metallic (TLL) phase is a many-boson state, espe- 
cially for small Cl, where the background heavily fluctu- 
ates. Note that this regime is not that easy accessible 
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